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Hannestad and Scherrer analysed self-interacting Warm and Hot Dark Matter and found power 
spectra in disagreement with previous papers. We argue that they include self-interactions in a weak- 
coupling approximation, as for photons after recombination. The "tightly coupled" approximation 
used for pre-recombination photons, and previous discussions of self-interacting HMD, would have 
been more appropriate. Their approximation generates a Boltzmann hierarchy leading to a stiff 
system of equations. Furthermore contamination by gauge modes could invalidate their results. 
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Neutrinos with strong self-interactions were originally 
suggested as a dark matter candidate]] by Raffelt and 
Silk ffl, who estimated that density perturbations in 
such a fluid would damp by diffusion ("Silk damping"), 
rather than by free-streaming. In a subsequent analysis — 
partially numerical, partially analytic — Atrio-Barandela 
and Davidson (ABD) g computed a power spectrum in 
agreement with this expectation. More recently, power 
spectra were computed by Hannestad and Scherrer || for 
self-interacting hot and warm dark matter, using CMB- 
FAST 0| . They find that the shape of the power spec- 
trum is essentially independent of the strength of the 
dark matter self-interactions, but depends on the mass 
of the dark matter particle, suggesting that the damp- 
ing of density fluctuations is the due to free-streaming 
rather than diffusion. Their power spectra do not show 
the damping due to diffusion and the short distance os- 
cillations that one would expect for an interacting dark 
matter fluid [0J2). We believe this is due to an inappro- 
priate analytic approximation, and possibly to contami- 
nation by gauge modes. 

HS study the evolution of perturbations ^>(x l , ft, r) in 
the dark matter phase space distribution /: 



£[/] = c[f] 



(3) 



/(x\ft,"r) = / (l + ft, T)) 
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/o is the homogeneous and isotropic equilibrium phase 
space distribution, and x 1 ,t are comoving coordinates 
(see [U for notation). In the absence of interactions, / 
satisfies the Boltzman equation for free particles 



or ot ox 1 or oqj 
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The CMBFAST Boltzman code solves this equation for 
'J to evolve the HDM density perturbations. 

The Boltzman equation for an interacting species is 



*a light bosonic particle could behave very differently Q. 



where the collision term C[f] is difficult to solve. 
Schematically: 

C[.fa] = 2^/ dn b dU c dU d \M\ 2 S(f a f b - f c f d ) (4) 

where dH b = d 3 p b /[(2ir) 3 2E b ], M. is the matrix element 
for the process a + b <-» c + d, 6 is the 4-momentum con- 
servation delta-function, and we have assumed low occu- 
pation numbers. An equilibrium distribution is defined 
to satisfy C[f] = 0, so C[fo] = = C[/o], and 



(5) 



as noted by HS. We first discuss an approximate solu- 
tion suitable for strongly interacting particles, and then 
compare to the method of HS. 

The timescale associated with the LHS of equation (|^) , 
for horizon-sized modes, is the age of the Universe tjj- 
The timescale associated with the RHS is the interaction 
timescale r„ and t v <C tjj for strong self-interactions. 
In this case, (||) can be solved in the "tightly-coupled" 
approximation used for pre-recombination photons: per- 
turbations will be very close to locally in equilibrium, 
because the interactions are fast. So 

f{x\ qi ,T)=gW{x\T)+g^{x\q u T) (6) 

where 

g^{x\T)=eM{E-v{x\T))/T{x\T)} , (7) 

is the locally in equilibrium distribution, so is the 
perturbation away from local equilibrium. By definition 
C^ -*] = 0, so the Boltzman equation is 



(8) 



This implies g^ ~ (tv/t^/o^. The equality in eqn 
(Q) is the analytic approximation required to generate 
a hierarchy of Boltzmann equations similar to the one 
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solved by HS. From the approximation in eqn flq), we 
showed that density perturbations in interacting HDM 
models would be damped by diffusion B, as anticipated 
in 0. (This is also the result obtained by HS in their 
Appendix 2.) 

HS used a different approach from us; they take 
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and generate a Boltzmann hierarchy. The resulting equa- 
tions are, however, very stiff. Their approximation (equa- 
tion |^) diverges in the limit of very large scattering cross 
sections and will not model the fact that fast interac- 
tions can bring ^ into a local equilibrium form almost 
satisfying C[^] = 0. Therefore, we believe that equation 
(0) does not correctly model the HDM self-interactions, 
and this could explain why HS's power spectra look like 
non-interacting Hot Dark Matter. Equation (0) would 
be valid in the weak coupling approximation used for 
photons after recombination, as opposed to the pre- 
recombination tight coupling approximation leading to 
equation (||). 

Note that both ABD and HS work in synchronous 
gauge, where density perturbations in the CDM regime 
are defined in the matter rest frame. However, the di- 
vergence of the fluid velocity is non-zero for relativistic 
neutrinos, so one must verify that U = after the v 
become non-relativistic, to ensure that unphysical gauge 
modes do not contaminate the solution. This could be 
problematic in the full numerical approach of HS. 

As mentioned, HS discuss the relation of their paper to 
ABD in an Appendix, where they reproduce the analytic 
estimates showing that the perturbations should os- 
cillate within the horizon and be damped on the diffusion 
scale. They do not comment on the discrepancy between 
the estimates and their plots. They suggest that the dif- 
ference between the ABD P(k) plot and theirs is due to 
an extra HO term appearing in our equations, where 6 
is the divergence of the fluid velocity. This term is due 
to a difference in conventions: as stated in our paper, 
we use the conventions of || for our analytic derivations, 
so 8 = ik ■ v/a, and our equations are correct (see eqn 
85.8 in j(|). HS use the conventions of ||, according to 
which = ik ■ v, and the offending HO term does not 
appear. In MB notation, where <5 represents the density 
perturbation, equations ABD [40] and [41] are: 



regime density perturbations grow like Cold Dark Mat- 
ter. However, HS are correct to doubt our plots; our 
treatment of the transition from relativistic neutrino fluid 
to non-relativistic fluid was simplistic, and we did not 
correctly take into account the effect of viscosity on small 
scales (a factor 3 was missing in our programme) . During 
the transition period from radiation to matter dominated 
regimes we assumed that w — P/ p, was linear with con- 
formal time. Other relationships gave very similar power 
spectra. We have corrected the viscosity term, and plot 
the power spectrum in figure [l]. This correction has a 
negligeable effect: the small scale oscillations are slightly 
reduced with respect to j^], and our results agree with 
expectations. 




FIG. 1. Power spectrum for two different "sticky neu- 
trino" models. The solid lines correspond (in decreasing am- 
plitude) to interactions with co-moving mean free paths at 
T 7 = 10 eV of .01 Mpc and .1 Mpc. These give rise to esti- 
mated damping scales of 7 Mpc -1 and 3 Mpc -1 , respectively. 
The dashed line correspond to standard CDM and is plot- 
ted for comparison. The y-axis scale is arbitrary. We took 
H = 50Km s^Mpc" 1 . 
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where w = P/p, P = pressure, p — density, a is the scale 
factor, a the viscosity, k the wavelength in conformal 
units, h and ry are the metric perturbations, and deriva- 
tives are with respect to conformal time. In the relativis- 
tic regime w = 1/3, w = and neutrinos behave like a 
damped harmonic oscillator, while in the non-relativistic 
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